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Abstract 

In this paper we examine Grosswald’s conjecture on g{p), the least primitive root mod¬ 
ulo p. Assuming the Generalized Riemann Hypothesis (GRH), and building on previous 
work by Gohen, Oliveira e Silva and Trudgian, we resolve Grosswald’s conjecture by 
showing that g{p) < ^fp — 2 for all p > 409. Our method also shows that under GRH 
we have g{p) < y/p — 2 for all p > 2791, where g{p) is the least prime primitive root 
modulo p. 


1 Introduction 

Let g{p) denote the least primitive root modulo p. Burgess [2] showed that g{j)) -C 
for any e > 0. Grosswald [1] conjectured that 

g{p) < y/p-2, 
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for all primes p > 409. Clearly, Burgess’ result implies ([T]) for all sufficiently large p. In [3] it 
was shown that ([T]) is true for all 409 < p < 2.5 • 10^® and for all p > 3.38 • 10^^. In this paper, 
contingent on the Generalized Riemann Hypothesis (GRH) we prove ([T]) for the remaining 
values of p. In fact, we prove a stronger result. Let g{p) denote the least prime primitive 
root modulo p. 


Theorem 1. Assume GRH. Then g{p) < y/p — 2 for all primes p > 2791 and g{p) < y/p — 2 
for all primes p > 409. 


We use Theorem [T] to make the following improvement to Theorem G in [1] . 


Corollary 1. Assume GRH. For all primes p, the modular group r(p) can he generated by 


the matrix 


1 p 
0 1 


and p{p — l)(p + 1)/12 canonically chosen hyperbolic elements. 


Proof. Our Theorem [T] covers the range of p in Theorem 2 in [4] , whence the result follows. □ 

The layout of this paper is as follows. In ^we prove an explicit bound on the least prime 
primitive root. Using this, we are able to prove Theorem [1] for all p > 10^^. We introduce a 
sieving inequality in ^ which allows us to complete the proof of Theorem [H 

Throughout the paper we write u:{m) to denote the number of distinct prime divisors 
of m. 


2 An explicit bound on the least prime primitive root 

Theorem 2. Assume GRH. When p > 10®, the least prime primitive root g{p) satisfies 


g{p) < 0 {‘2‘^^ - l) logp^ 


( 2 ) 


Before we prove Theorem [2] we remark that the bound in ([2]) is not the sharpest known. 
Shoup [7] has proved that 

g{p) < uj{p - l)^(log(a;(p - 1)) + l)^(logp)^ 


and, as remarked by Martin |5l p. 279], this bound also holds for g{p). While this su¬ 
persedes (j2]) for all sufficiently large primes p, the utility of Theorem [2] is in providing a 
completely explicit estimate. 

We require the following result, which is easily established following the proof of Lemma 2.1 

in d]. 


Lemma 1 (Bach). Let x(n) denote a non-principal Dirichlet character modulo p. 
X >1, we have: 


n<x 


< ^ + 3 , 
- 20 


When 

( 3 ) 
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and 


where 


^A(n)x(n) (l - -) 

\ X' 

n<x 


< c{p,x)y/x\ogp, 


c{p, x) 


2 / logx + 2 

3 V ^/x x^/'^ J V \ogp) ^/x\ogp 


( 4 ) 


Proof of Theoreml^ We may assume g{p) > 1099, or else there is nothing to prove; indeed, 
one has ((8/5) log(10®))^ ~ 1099.4. As in Burgess [21 §4], we rewrite the function 


fin) = 


1 if n is a primitive root modulo p 
0 otherwise. 


as 


fin) = 


0(P - 1) 

p — 1 


d\p-l ^ X 
d>l 


n) 


where 'Yfx f^iken over all Dirichlet characters modulo p of order d. Suppose that f{n) = 0 
for all primes n (and hence all prime powers) with n < x. We may assume x > 1099. We 
have 

(5) 

d\p-l ' X 
d>l 

for all prime powers n < x. We multiply ([5]) by A(n)(l — n/x) and sum over all n < x. This 
gives 

/i(d) 


n<x 


5:Aw(i-/) + Eg|EEAWxW(i-;)^o, 

d\p—l ^ ^ ^ X 
d>l 


We apply Lemma [Hand observe that c{p, x) < 7/9 provided p > 10® and x > 1099. Therefore 


d\p-l 

d>l 

tJ.{d)jtO 


which implies 

- 1) logp < 5 - 1) logp- 


The resnlt follows. □ 


Corollary 2. Theorem [I] is true except possibly when p G (2.5 • 10^^, 10^^) and ujfp — 1) G 
[7,28]. 
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Po 

102 

10^ 

10® 

10® 

lO'o 

C 

2.1127 

1.6821 

1.5556 

1.496 

1.4614 

Po 

10^2 

10'" 

10'® 

0 

00 

1020 

C 

1.4389 

1.4231 

1.4114 

1.4023 

1.3952 


Table 1: Values of C for various choices of po 


Proof. Unconditionally, Robin (see |6l Thm 11]) proved 

logn 


uj{n) < 1.385 


log log n ’ 


(n > 3). 


This, when combined with Theorem [2l shows that g{p) < y/p — 2 provided p > 10^®. Hence 
we may assume p < 10'^®. But the assumption p < 10^® leads to u:{p — 1) <31 and 
now Theorem [2] implies the result provided p < 10^^. Repeating this process, we find that 
u:{p — 1) < 28 and p < 10^^. On the other hand, in light of the computations carried out in 
0 §4|, we may assume p > 2.5 • 10^®, in which case Theorem [2] proves the result provided 
(^{p — 1 ) < 6 . □ 


3 Computations using the sieve 

In this section we follow closely the argument given in [3l §3]. Suppose e is an even divisor of 
p — 1. Let pi,... ,Ps he the primes dividing p — 1 that do not divide e. Set 5 = 1 — 
and set n = u:{p — 1). In applying our method it is essential to choose 5 > 0. 

Theorem 3. Assume GRH. If g{p) > x, then we have: 

We postpone the proof of Theorem [3] until 1 13.21 From Theorem [3] we immediately obtain 
the following corollary which is more readily applied. 

Corollary 3. Assume GRH. If p> po, then 

Hp) <R ( 2 + ^^'\ 2"-‘ log pi , (7) 

where the constant is given in TableUi 


g{p)< 2c{p,x) 



Proof. When p > pg, the right-hand side of ([7]) is bounded below by x := ((4/3)4 logpo)^- 
Hence we may assume g(p) > x, or else there is nothing to prove. Now Theorem [3] establishes 
the result with C(po) := 2c(po,x). □ 
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Our proof of Theorem [T] will apply Theorem [3] directly, but we have included Corollary [3] 
as it may have application elsewhere. 

Corollary 4. Theorem [I] is true except possibly when uj{p — 1) G [12,13,14], 

Proof. In light of Corollary |2l we may assume 7 < n < 28 and 2.5 ■ 10^® < p < 10^^. We 
have the obvious lower bound p — 1 > qi.. .qn, where g* denotes the ith prime, and hence 
we may assume 

p > max i 1 + g*, 2.5 ■ 10^® 

I i=i 

For example, when n = 15, this leads to p > 6.1 ■ 10^^. 

Set X = — 2. We may assume g(p) > x or else there is nothing to prove. Hence the 

conclusion of Theorem [3] holds. Now each choice of s allows us to show that 

g{p) < y/p-2 (8) 

holds when p is larger than an explicitly computable value; one simply bounds the right-hand 
side of ([6]) from above, usin^ h > 1 — Yl^=n-s+iPi- course, we then choose the value of s 
that gives the best result. For example, when n = 15 we hnd that s = 12 leads to 

5 > 1 - (1/7 + 1/13 + ■ ■ • + 1/47) > 0.3717 

and therefore the right-hand side of (j6]) is less than — 2 provided p > 3.2 ■ 10^®; hence 
any exception to (IH]) must satisfy p < 3.2 ■ 10^®. Notice that because our lower and upper 
bounds on any potential exceptions overlap, this proves the result when n = 15. The best 
choice turns out to be s = n — 2 when 7 < n < 12 and s = n — 3 when 13 < n < 28. In fact, 
this is enough to prove ([8]) except when n = 12,13,14. The lower bound of 1 -|- nr=i 
the job when 15 < n < 28 and the lower bound of 1.6 • 10^® does the job when 7 < n < 11. 

□ 



3.1 An algorithm 

In order to deal with the cases when n = 12,13,14, we introduce an algorithm. Before giving 
the algorithm, we explain the main idea. 

Suppose n = 14. Using the idea presented in the proof of Corollary 01 we find that 
any exception to Grosswald’s conjecture must lie in the interval (1.30 ■ 10^®, 1.71 ■ 10^®). In 
principle one could check the conjecture directly for each prime p in this interval, but the size 
of the interval makes this prohibitive. There are 2.05 ■ 10^® odd values of p to consider. (Of 
course many of these are not prime, but there are still about 10^“^ primes in this interval.) 
Instead, we break the problem into cases depending upon which primes divide p — 1. We 
already know that 2 divides p — 1, so we start with the prime 3. If 3 divides p — 1, then 

^ One helpful fact — if the right-hand side of ([S]) is less than — 2 for some p, then the same is true for 
all larger p. 
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we have one third as many values of p to check, roughly 6.83 • values of p. If 3 does 
not divide p — 1, then this leads to an improved lower bound on p, as well as an improved 
lower bound on 6 and hence an improved upper bound on p; in short, the interval under 
consideration shrinks. In this particular case, the interval shrinks to (2.04 ■ 10^^, 1.45 • 10^®), 
which is empty, so there is nothing to check. 

Returning to the case where 3 divides p — 1, the number of exceptions is still quite large. 
However, we may consider whether 5 divides p — 1. We continue in this way until the number 
of possible values of p under consideration is less than 10®. At that point we go through the 
list and throw out all values of p except those where p is prime, a;(p — 1) = 14, and p — 1 
satisfies the given divisibility conditions (depending upon which sub-case we are considering). 
We append these exceptional values of p to a list and continue this recursive procedure until 
we have exhausted all possibilities. One can easily find the least prime primitive root for the 
list of exceptions via standard methods and check the conjecture directly. When n = 14, this 
algorithm takes only 7 seconds (on a 2.7 GHz iMac) to complete and the list of exceptions 
is empty, so there is nothing further to check. The number of exceptions for other values of 
n is given in Table |2l 


Table 2: Number of exceptions for n = 12,13,14 


n 

12 

13 

14 

of exceptions 

61,114 

6,916 

0 


For completeness, we give the pseudocode for our recursive algorithm. Suppose X UY 
consists of the first k primes for some k G Z>o. Algorithm [1] will verify g{p) < y/p — 2 
when a;(p — 1) = n under the assumption that q divides p — 1 for all g G X and q does not 
divide p — 1 for all g G H. More precisely, rather than verifying the conjecture for all p, 
the algorithm will generate a manageable list of possible exceptions which can be checked 
individually, as described above. The sets X and Y are allowed to be empty, although in 
practice we may always assume 2 G X. (Running the algorithm with n = 14, X = {2}, 
X = 0 will carry out the computation described above.) 

Proof of Theorem [H We have implemented Algorithm [T] in Sage. Running our code on 
n = 12,13,14, including finding the least prime primitive root and checking the conjecture 
directly for the list of 68,030 exceptions, takes about 4.5 minutes. In light of Corollary 01 
this proves the theorem. □ 


3.2 Proof of Theorem [3] 

Let p be an odd prime. Let e be an even divisor of p — 1. We say that n is e-free if the 
equation = n (mod p) is insoluble for all divisors d of e with d > 1. An integer is a 
primitive root if and only if it is (p — l)-free. We define the function 


fe{n) 


1 if n is e-free 
0 otherwise. 
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Algorithm 1 Grosswald(n,X,Y) 

1: L ■= first n primes not in Y 

2: lower := max{pro(iMct(L) + 1, 2-10^^} 

3: upper := 0 

4: for s G {1,..., n — 1} do 
5: M := largest s primes in L 

5 := 1 — XlgeM q 

7: if 5 < 0 then 

8: continue 

9: end if 

10: Choose p large enough so that when x = ^ — 2, we have: 

(^2c{p, x) (^2 + log(p)^ < x 

11: if upper = 0 or p < upper then 

12: upper := p; s := s; (5 := S] M := M 

13: end if 

14: end for 

15: print(n, s, S, M, lower, upper) 

16: prodX := product (X) 

17: enum := {upper — lower)/prodX 

18: if enum < 0 then 

19: print (“Nothing more to check.”) 

20: else if enum > 10® then 
21: print(“Break into cases.”) 

22: q :=smallest prime not in X U Y 

23: Grosswald{n, X U {g}, Y) 

24: Grosswald{n, X,Y U {g}) 

25: else 

26: print (“Find all exceptions.”) 

27: / := [{lower — 1) /prodX, {upper — 1) /prodX] flZ 

28: Sieve out elements of I divisible by primes in Y 

29: for /c G / do 

30: p := k* prodX + 1 

31: if is_pseudoprime(p) and length(prime_divisors(p — l))=n then 

32: append p to global list of exceptions 

33: end if 

34: end for 

35: end if 
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Define the multiplicative function 6{n) = 0(n)/n, where 4>{n) is Euler’s totient function. We 
rewrite feiji) as 


We see 


Thus 


d>\ 


'^fpAn) - (s-l)/e(n) 

i=l 


= 1 if n is (pje)-free for all i 
< 0 otherwise. 


fp-i{n) > ^ /p,e(n) - (s - l)/e(n) 

i=l 
s 

= X] (fpie{n) - 9{pi)fe{n)) + 6fe{n). 


Observe that 


2=1 


fpAn) - 0{pi)fe{n) = e{pie) ^ ^ ^ x{n) 

d\pie ^ Xd 
d\e 
d>l 




4^ 


Inserting (ITOll into ([9]) leads to 


1 + i E OM E 7r§ E x(")+E E x(n) < 

d>l 


2=1 


/p-iW 

59{e) ' 


(9) 


( 10 ) 


( 11 ) 


Suppose fp-i{n) = 0 for all primes n (and hence all prime powers) with n < x. We multiply 
(ITTl) by A(n)(l — n/x), sum over all n < x, which yields 


Eaw(i-9 


n<x 


E^wE^EEaWawIi- 


5 


2=1 




E^EExWA(")(i-;)<0. 

d\e ^ ^ X 
d>l 










We write 


5:A(n)(l-^)=2 + GW. 

n<x 

and use the estimates in ([3]) and ([1]) to obtain 


f + G{x) 



S 

c(p, x)^/x\ogp + 

d\e 

fi{pid)y^0 


c{p,x)y/xlogp 


d\e 

d>l 

tJ.{d)^0 


< c{p, x)^logp 


= c{p, x)y/x\ogp 



This leads to 


X < 


2|G(x)| 


X 


+ 2c(p, x) 


16 , , 
< T7T + ^ + 2cP,x 

10 wx 


2 + -—^ 


2 + f-^ 1 nn-s 


logp 

logp 


< 2c(p,x)l2 + ^]2^-^ 


log p -^-log p 

90 3 


The result follows. 


□ 
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